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Introduction 

The  Chemical  Oxygen  Iodine  Laser  (COIL)  is  the  only  demonstrated  example  of  a 
chemically  pumped  electronic  transition  laser ^  COIL  systems,  operating  at  a  wavelength 
of  1.315  p,m,  provide  very  high  energy  laser  beams.  These  devices  normally  function  in  a 
continuous  mode,  but  can  also  be  modulated  rapidly  by  Zeeman  gain-switching2-4. 
Recently,  frequency  doubling  of  a  COIL  device  has  been  achieved^.  Several  military  and 
commercial  applications  for  COIL  systems  are  being  considered.  Experimental  programs 
aimed  at  increasing  the  efficiency  of  the  oxygen-iodine  laser  and  raising  the  working 
pressure  of  the  device  are  currently  in  progress.  In  parallel  with  these  efforts, 
computational  models  are  being  developed  so  that  energy  loss  mechanisms  may  be 
identified,  and  methods  for  eliminating  or  reducing  these  problems  may  be  explored^.  Any 
realistic  model  of  the  laser  can  be  divided  into  physical  and  chemical  components.  The 
physical  part  of  the  problem  deals  with  reagent  mixing  and  transport  properties,  while  the 
chemical  component  deals  with  reactions  that  generate,  excite,  and  deactivate  atomic 
iodine^  16.  At  present,  the  reliability  of  computer  modeling  is  compromised  by  unresolved 
questions  regarding  the  I2  dissociation  mechanism,  and  the  lack  of  rate  constant  data  for 
several  key  processes. 

In  COIL  systems,  I2  is  dissociated  by  multiple  collisions  with  electronically 
metastable  species.  It  is  thought  that  the  two  most  important  metastables  in  this  process  are 
02(alA)  and  I(^Pi/2).  However,  the  mechanism  by  which  I2  is  dissociated  in  l2/02(alA) 
mixtures  has  been  the  subject  of  several  investigations,  and  of  some  controversy^'l^. 
Despite  these  efforts,  the  mechanism  is  still  not  fully  understood.  As  02(alA)  does  not 
posses  sufficient  energy  to  dissociate  I2  in  a  single  collision,  sequences  involving  energy 
pooling  and  multiple  energy  transfer  events  have  been  considered.  Current  theories^'^^ 
suggest  that  the  first  step  may  be  the  electronic  to  vibrational  transfer  (E-V)  process 

I2(X)  +  02(alA)  ^  I2(X)  v"=33-44  -h  02(X)  (1) 

Subsequent  collisions  of  02(a^A)  with  I2(X)  v">20  can  result  in  either  direct  dissociation, 
or  electronic  excitation  of  I2  to  the  A'(2u),  A(lu),  or  BCOu"^)  states.  There  are  then  several 


reactions,  generalized  by  the  equations 

02(alA)  +  l2t->  21  +  02(X)  (2a) 

02(alA)  + 12*^  21  +  02(X)  (2b) 

that  can  contribute  to  the  dissociation.  In  these  equations  12'!'  denotes  I2(X)  v">20,  and  I2* 
is  electronically  excited  I2.  Once  atomic  I  has  been  liberated,  it  is  rapidly  excited  by  the 
process 


.on/",  i 


Avallabilitf  Cod^S 
Avail  aad/er 
Lst  Spsetal 


02(alA)  +  I(2P3y2)^  I(2Pl/2)  +  02(X)  (3) 

In  addition  to  being  the  species  responsible  for  lasing,  there  are  indications  that  I(^Pi/2) 
also  participates  in  the  dissociation  via  the  reactions 

I(2Pl/2)  +  l2^  l2^  +  I(2P3/2)  (4) 

and,  possibly  I(^Pl/2)  +  ^2^  31  (5) 

The  most  extensive  investigation  of  the  dissociation  mechanism  was  conducted  by 
Heidner  et  al7,  who  used  chemiluminescence  to  follow  the  concentrations  of  I(^Pi/2). 
12(A),  12(B),  02(alA),  and  02(bl5:‘^)  in  a  discharge  flow-tube.  This  was  the  first  study  to 
show  that  slow  "initiation"  reactions  preceded  rapid  I2  dissociation,  and  that  an  unidentified 
excited  state  of  I2  was  involved  in  the  dissociation.  Two  distinct  kinetic  models  were 
developed  and  evaluated  by  Heidner  et  al.7  These  differed  considerably  in  the  rate 
constants  they  assumed  for  generation  and  removal  of  excited  I2.  Neither  model  gave  an 
entirely  satisfactory  fit  to  the  experimental  data.  Heidner  et  al.^  considered  the  possibilities 
of  the  unidentified  excited  species  being  either  l2(A')  or  12"^.  Problems  were  noted  for 
either  assignment.  The  kinetics  indicated  that  the  excited  species  was  fairly  stable,  and 
almost  immune  to  quenching  by  Ar.  On  the  basis  of  physical  intuition  regarding  vibrational 
relaxation,  Heidner  et  al.7  felt  that  the  stability  of  the  excited  species  was  more  consistent 
with  I2(A').  They  were  uncertain  of  this  assignment,  however,  because  the  A’  state  lies  too 
far  above  thermally  populated  I2(X)  levels  (Te(A’)=10042  cm*^)  to  be  excited  by  a  single 
collision  with  C)2(alA). 

The  stability  of  l2^  was  probed  in  subsequent  experiments.  Hall  et  al.l®  reported  a 
relatively  large  rate  constant  of  2.3xl0'll  cm3  g-l  for  deactivation  of  l2(X,  v"=40)  by  Ar. 
This  result  was  not  necessarily  considered  to  be  in  conflict  with  the  observations  of 
Heidner  et  al.^  Reactions  1  and  4  populate  vibrational  levels  around  v"=40,  whereas 
transfer  from  02(alA)  or  I(2Pi/2)  can  dissociate  12"*^  from  levels  above  v"=20.  Therefore, 
many  collisions  may  be  needed  to  relax  the  nascent  l2^  to  levels  below  v"=20.  Indeed, 
David  et  al.^  have  shown  that  the  flow-tube  kinetics  can  be  modeled  using  the  rate 
constants  of  Hall  et  al.^O  Constraints  required  by  the  model  were  that  relaxation  proceeds 
via  sequential  v— >v-l  steps,  and  that  the  rate  constants  for  these  steps  scale  linearly  with  v. 

The  models  proposed  by  Heidner  et  al.^  were  intended  to  stimulate  further  research. 
They  did  not  fully  reproduce  the  experimental  data,  and  were  not  unique.  Even  so,  these 
models  have  provided  estimates  for  several  rate  constants  that  have  not  been  examined  by 
other  means.  Consequently,  computer  models  of  COIL  systems  that  have  been  developed 
at  the  Air  Force  Phillips  laboratory,  and  by  other  research  groups^,  have  relied  heavily  on 
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the  flow-tube  rate  constants.  No  doubt,  some  the  shortcomings  of  the  laser  models  can  be 
traced  to  inaccuracies  in  the  rate  constant  data  base.  Research  carried  out  under  grant 
AFOSR  49620-92- J-0073  was  primarily  directed  towards  the  measurement  of  energy 
transfer  processes  of  relevance  to  COIL  systems 

1.  Spectroscopy  and  Electronic  Energy  Transfer  Dynamics  of  I2. 

A.  Spectroscopic  studies  of  the  I2  D-X  transition. 

To  perform  energy  transfer  measurements  on  l2(X)  v">20  we  needed  a  sensitive, 
state-specific  method  for  observing  these  highly  excited  vibrational  levels.  Laser  induced 
fluorescence  (LIF)  of  the  intense  D-X  transition  of  I2  proved  to  be  ideal  for  this  purpose, 
but  it  quickly  became  apparent  that  extension  of  the  spectroscopic  data  base  was  an 
essential  prerequisite  for  energy  transfer  studies.  Problems  traceable  to  the  spectroscopic 
constants  were  encountered  by  Nota  et  al.^^  when  they  tried  to  extract  population 
distributions  fi'om  D<— X  spectra  recorded  in  02(alA)/l2  mixtures.  The  available  constants 
were  not  accurate  enough  for  direct  inversion  of  the  line  intensity  data.  Due  to  large 
differences  between  the  X  and  D  state  equilibrium  intemuclear  distances,  transitions 
suitable  for  observing  I2(X)  levels  in  the  range  20<v"<50  terminate  on  the  0^v"<40 
vibrational  levels  of  12(D).  Prior  to  the  present  studies,  the  most  accurate  data  for  the  lower 
levels  of  12(D)  had  been  obtained  from  optical-optical  double  resonance  (OODR) 
measurements^^.  The  bands  of  interest  for  probing  l2^  were  not  well-characterized  by  this 
approach. 

Two  techniques  were  used  to  examine  the  12(D)  (4<v'<37)<— X(23<v"^2)  bands. 
In  one  series  of  experiments,  photofi'agmentation  of  l2Arn  clusters  was  used  to  populate 
excited  vibrational  levels  of  the  ground  state^^  (^2^  was  ejected  from  the  fragmenting 
clusters,  see  below).  Very  low  rotational  temperatures  were  attained,  yielding  relatively 
uncongested  spectra.  In  the  second  series  of  experiments,  two  pulsed  lasers  populated  a 
specific  ro-vibrational  level  (vi  Ji)  of  I2(X)  by  stimulated  emission  pumping  (SEP).  A  third 
laser  was  then  used  to  observe  D-X  lines  originating  from  vj  Jj,  and  nearby  levels  that  were 
populated  by  collisions.  This  experimental  configuration  (shown  in  Figure  1)  was 
developed  for  studies  of  inelastic  collisions,  but  it  also  provided  a  convenient  method  for 
completing  the  spectroscopic  ground-work.  In  all  of  the  measurements,  D-X  spectra  were 
excited  by  the  frequency  doubled  output  from  a  tunable  dye  laser.  Absolute  wavenumber 
calibration  was  achieved  by  using  the  fundamental  of  the  dye  laser  to  simultaneously  record 
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Figure  1.  Experimental  apparatus  for  studies  of  I2  spectroscopy  and  energy  transfer 


the  I2  B-X  spectrum.  The  positions  of  the  B-X  lines  were  taken  from  the  work  of 
Gerstenkom  and  Luc^^. 

Rotationally  resolved  spectra  were  recorded  for  the  D-X  bands  12>v'>5  <-  v"=38 
and  37>v'>30  <—  v"=23.  Band  origins  and  rotational  constants  were  obtained  by  fitting 

the  line  positions  to  the  expression 

V  =  V,, +  b;  J  (J  +1)  -  B';r  (J"+l)  (1) 

(trial  fits  to  an  expression  that  included  the  centrifugal  distortion  constants  did  not  yield 
statistically  significant  values  for  these  parameters).  The  lower  state  rotational  constants 
were  found  to  be  in  good  agreement  with  those  obtained  by  Bacis  et  al.^®  and  Tellinghuisen 
et  al.21  As  these  earlier  studies  of  the  ground  state  were  of  greater  accuracy  than  our 
measurements,  the  final  fits  to  equation  (1)  were  made  with  the  lower  rotational  constants 
fixed  at  the  values  given  by  Bacis  et  al.^®  This  scheme  yielded  D  state  rotational  constants 
that  were  in  excellent  agreement  with  the  values  reported  by  Ishiwata  and  Tanaka^ 
Vibrational  term  energies  for  the  D  state  were  calculated  by  adding  the  ground  state 
vibrational  energies^O  to  the  fitted  band  origins.  The  resulting  vibronic  term  energies 
(relative  to  the  ground  state  minimum)  are  listed  in  Table  1.  The  new  term  energies  were 
compared  with  energies  predicted  using  the  constants  of  Ishiwata  and  Tanaka^7  At  v"=34 
the  predicted  energies  were  in  error  by  0.36  cm*l.  This  error  increased  with  decreasing  v", 
reaching  a  value  of  1.77  cm'^  at  v'=5.  Comparisons  with  energies  predicted  from  the 
constants  of  Bartels  et  al.22,  and  Hoy  and  Lipson^^  revealed  essentially  the  same 
discrepancies.  To  put  this  in  perspective,  a  predictive  accuracy  of  ±0.08  cm'^  is  needed  for 
the  inversion  of  line  intensity  data  from  congested  spectra. 

To  obtain  improved  spectroscopic  constants  for  12(D),  we  combined  our  results 
with  the  best  available  data  for  the  levels  above  v'=37.  Ishiwata  and  Tanaka^  noted  that 
their  data  were  most  accurate  for  the  range  76<v'<124.  For  these  levels  they  reported  that 
the  standard  deviation  of  a  global  fit  was  0.052  cm'l.  Vibronic  term  energies  for  the  ranges 
125<v'^01  and  202<v'^90  were  taken  from  the  studies  of  Bartels  et  al.22  ^nd  Hoy  and 
Lipson23,  respectively.  The  complete  set  of  vibronic  term  energies  were  fit  to  the  standard 
Dunham  expression 

T'’(v)  =  2;y,,„(v  +  1/2)” 

n 

The  data  were  weighted  according  to  the  inverse  squares  of  their  standard  deviations 
(1/(0. 15  cm-l)2  for  v'<124;  l/(0.5cm-l)2  for  v'>125).  A  satisfactory  fit  was  achieved 
using  a  ninth  order  polynomial.  As  the  constants  obtained  were  highly  correlated,  they 


4 


Table  1 


Vibrational  Term  Energies  for  12(D) 


v' 

v(cm-l) 

5(cm-l> 

4 

41451.97 

0.15 

5 

41545.84 

0.13 

6 

41639.56 

0.13 

7 

41733.01 

0.14 

8 

41826.28 

0.08 

9 

41919.19 

0.13 

10 

42011.95 

0.15 

11 

42104.47 

0.14 

30 

43810.37 

0.15 

31 

43897.37 

0.14 

32 

43984.05 

0.11 

33 

44070.47 

0.09 

34 

44156.55 

0.09 

35 

44242.35 

0.11 

36 

44327.85 

0.14 

37 

44413.09 

0.13 

♦ 

a.  1-0  errors  from  fits  to  equation  1. 


were  subjected  to  a  systematic  rounding  procedure  that  preserves  their  accuracy^'*.  The 
final  constants  are  given  in  Table  2,  where  they  are  compared  with  several  previous 
estimates.  The  residuals  of  the  fit,  determined  for  the  rounded  constants,  are  plotted  in 
Figure  2.  The  1-g  statistical  error  band  for  the  fitted  T^(v)  function  is  shown  in  Figure  3. 
Differences  between  the  D  state  energies  calculated  from  the  present  results,  and  those  from 
the  constants  of  Ishiwata  and  Tanaka^^  (5TD(v)=Tl^(v)present-T^(v)ref.  i?)  are  also  plotted 
in  Figure  3.  Local  fitting  problems  are  primarily  responsible  for  the  deviations.  The 
previous  constants  work  well  over  the  ranges  of  levels  ft’om  which  they  were  defined.  For 
the  purpose  of  extracting  vibrational  populations  from  D-X  band  intensity  data,  we  have 
calculated  a  new  RKR  curve  for  the  D  state,  and  a  matrix  of  the  relevant  D-X  Franck- 
Condon  factors.  These  results  are  listed  in  the  appendix  of  this  report. 

B.  Spectroscopy  and  relaxation  of  l2Ar„  clusters.  Geminate  recombination  and  cluster 
fragmentation. 

l2Arn  clusters  were  studied  in  a  free-jet  expansion.  When  excited  to  states  that 
correlated  with  the  first  ion-pair  manifold  of  I2,  the  clusters  yielded  visible  and  near  UV 
emissions.  The  ion-pair  states  were  accessed  by  193  nm  excitation  of  ground  state 
clusters,  or  near  UV  excitation  of  electronically  metastable  l2(A'  or  AjAr^.  Comparisons 
of  the  cluster  and  12/Ar  matrix  spectra  suggested  that  such  excitations  result  in 
fragmentation  and  "melting"  of  the  clusters  prior  to  emission. 

532  nm  photodissociation  of  I2  within  the  clusters  was  followed  by  geminate 
recombination  that  populated  the  A',  A,  and  X  states.  The  probability  for  ejection  of  the 
recombined  I2  from  the  clusters  during  the  relaxation  process  appeared  to  be  size 
dependent.  These  results  represent  tentative  steps  towards  the  study  of  I(^P3/2)  +  I(^P3/2) 
recombination  in  Ar  clusters.  They  also  indicate  that  dissociation,  recombination, 
fragmentation  processes  may  be  used  as  a  general  method  for  generating  metastable  species 
in  supersonic  expansions.  This  work  is  described  in  reference  18. 

C.  Energy  tranter  from  I(^Pii2)  to  hlX) 

Pulsed  photolysis  and  laser  probing  techniques  were  used  to  study  the  nascent  state 
distribution  resulting  fi-om  l2(X)  +  I(2Pi/2)  ->  l2^  +  I(^P3/2)-  Similar  experiments  had 
been  conducted  by  Hall  et  al.^®,  but  the  spectra  obtained  were  not  rotationally  resolved.  In 
our  recent  work,  l(^Pip)  atoms  were  generated  by  using  490  nm  light  to  photodissociate  I2 
via  the  B  state  continuum.  After  a  short  delay  (typically  l|is  for  a  100  mTorr  sample),  l2^ 
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Table  2 

Electronic  Term  Energies  and  Vibrational  Constants  for  12(D)  (cm-1) 
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a.  The  lo  uncertainties  for  these  constants  are  given  in  parentheses 

b.  Due  to  a  typographical  error  this  is  given  as  a  positive  constant  in  reference  17. 


Figure  3 


levels  populated  by  E-V  transfer  were  probed  by  recording  rotationally  resolved  D-X 
spectra. 

The  new  molecular  constants  for  the  D  state  were  used  to  determine  the  product 
state  distribution  from  an  inversion  of  the  Une  intensity  data.  Within  each  vibrational  level 
the  rotational  populations  were  well  represented  by  a  Boltzman  distribution  with  a 
temperature  of  310±20  K.  The  nascent  vibrational  state  population  distribution  is  shown  in 
Figure  4.  Note  that  the  distribution  peaks  at  v"=40,  which  coincides  with  near  resonant 
energy  transfer,  and  cuts-off  abmptly  above  v"=47.  The  distribution  obtained  from  these 
measurements  has  been  used  in  the  determination  of  l2^  deactivation  rate  constants  from 
state-to-state  transfer  rate  constants  (see  below). 

D.  Energy  tranter  from  electronically  excited  oxygen  to  l2(X) 

An  attempt  to  examine  l2(X)  +  02(al A)  -»  l2'^  +  02(X)  E-V  energy  transfer  was 
made.  For  this  purpose,  02(a^A)  was  generated  by  using  248  nm  pulses  from  a  KrF  laser 
to  photodissociate  O3  in  flowing  O3/I2  mixtures.  As  in  the  experiments  described  above,  a 
tunable  probe  laser  was  used  to  examine  the  I2*  product  state  distribution  by  exciting  the 
valence  to  ion-pair  transitions.  Unexpectedly,  the  probe  laser  spectra  were  dominated  by 
the  I2  D'-A'  bands.  Observation  of  these  bands  depended  on  the  photolysis  of  O3,  and  it 
was  verified  that  the  A'  state  was  not  being  populated  by  an  indirect  mechanism  involving 
248  nm  excitation  of  I2.  The  rate  at  which  the  A'  state  was  excited  was  consistent  with  a 
single  collision  event.  High  yields  of  0(^D)  atom  are  produced  by  photolysis  of  O3  at  248 
nm.  The  transfer  process  0(lD)+l2(X)-»  l2(A’)+0(3p)  could  account  for  the  appearance 
of  l2(A'),  but  it  is  unlikely  that  this  channel  is  significant.  It  is  known  that  O  atoms  react 
rapidly  with  I2  to  form  lO+I.  If  a  parallel  energy  transfer  channel  does  exist,  it  is  more 
likely  to  proceed  by  the  near-resonant  0(^D)+  l2(X)-»l2(B)-i-0(3p)  path.  An  alternative  A' 
state  excitation  mechanism  becomes  apparent  when  the  excess  energy  delivered  by  the 
photolysis  pulse  is  considered.  The  energetic  threshold  for  dissociation  of  O3  into 
0(lD)-i-02(alA)  is  32,260  cm"^  so  248  nm  light  provides  about  8062  cm'^  of  excess 
energy.  When  this  is  partitioned  between  the  recoiling  fragments,  the  02(alA)  fragment  will 
receive  about  2700  cm'^.  Combined  with  the  electronic  energy,  this  sufficient  to  open  the 
channel  02(alA)*+l2(X)-^l2(A’)+02(X). 
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Vibrational  level 


Figure  4.  Vibrational  population  distribution  resulting 


2.  Rotational  and  Vibrational  Relaxation  of  12'^ 


A.  Experiment  and  results 

Conceptually,  the  technique  used  to  study  inelastic  collisions  involving  I2'*’  was 
straightforward^^.  Figure  1  shows  one  variant  of  the  apparatus  used.  Pulsed  lasers  were 
used  to  excite  a  single  ro-vibrational  level  of  l2(X)  by  stimulated  emission  pumping  (SEP). 
After  a  suitable  delay,  molecules  in  the  initially  excited  level,  and  in  levels  populated  by 
collisions,  were  detected  by  a  third  laser  that  excited  the  D-X  transition.  A  small 
monochromator  (0.25  m)  was  used  to  selectively  detect  the  D-X  fluorescence. 

In  order  to  study  rotational  and  vibrational  energy  transfer,  a  small  amount  of  I2 
(typically  10-25  mTorr)  was  mixed  with  0.05-1  Torr  of  a  collision  parmer.  Probe  laser 
spectra  were  then  recorded  for  a  number  of  fixed  delays  between  the  state  preparation  and 
probe  laser  pulses.  The  energy  transfer  processes  examined  involved  I2  in  one  of  the  initial 
levels  (v"=23,  J"=57),  (v"=38,  J"=57),  or  (v"=42,  J"=17),  colliding  with  He,  Ar,  I2,  O2, 
N2,  CI2  or  H2O.  Self  transfer  was  examined  using  an  I2  pressure  of  200  mTorr.  Figure  5 
shows  a  typical  probe  laser  spectrum.  This  trace  was  recorded  for  a  mixture  of  20  mTorr 
of  I2  with  60  mTorr  of  H2O.  The  delay  between  the  dump  pulse  and  probe  pulses  was 
200  ns.  Lines  originating  from  levels  populated  by  rotational  energy  transfer  are  prominent 
in  this  trace.  Weaker  features  produced  by  v"=23— >22  relaxation  can  also  be  seen  (note 
that  vibrational  transfer  occurred  with  very  little  change  in  the  angular  momentum). 

For  each  combination  of  an  initial  level  and  a  collision  parmer,  the  total  removal  rate 
constant  was  determined  by  fixing  the  probe  laser  to  monitor  the  initial  level,  and 
monitoring  the  fluorescence  intensity  as  a  function  of  dump-probe  pulse  delay  time.  Good 
single  exponential  decays  were  observed.  Decay  rates  were  determined  by  fitting  to  the 
expression  I(t)  =  Ae'*^^  where  A  is  a  constant  and  k  =  ki2'[l2]  +  k'M[M].  ki2'  and  k'M  are 

the  total  transfer  rate  constants  for  collisions  with  I2  and  M,  respectively.  Rate  constants, 
derived  from  Stem-Volmer  plots,  are  listed  in  Table  3.  All  M^l2  rate  constants  were 
corrected  for  self-transfer. 

Rate  constants  for  state-to-state  ro-vibrational  energy  transfer  (vj  Jj— >vf  Jf)  were 
determined  from  rotationally  resolved  spectra  taken  under  single  collision  conditions.  An 
indirect  method  was  used  to  verify  that  multiple  collision  events  were  not  modifying  the 
intensity  distributions  seen  in  these  experiments.  Spectra  were  taken  at  a  variety  of 
pressures  and  delay  times  for  each  collision  partner,  and  conditions  where  the  intensities  of 
the  collisionally  populated  spectral  features  began  to  show  non-linear  time  or  pressure 
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2535.5  2536.0  2536.5  2537.0 

Wavelength  (A) 

Figure  5.  I2  D-X  spectrum  showing  levels  populated  by  rotational  and 
vibrational  energy  transfer 


Table  3 


Rate  Constants  for  Population  Loss  from  Individual  Ro-Vibrational  Levels 

0fl2(X) 


Collision 

partner 

kT(v=23,  J=57)a 

kT(v=38,  J=49)a 

kT(v=42,  J= 

He 

6.4 

5.4 

6.6 

Ar 

5.7 

5.5 

5.0 

l2 

8.7 

6.6 

8.6 

O2 

7.1 

7.0 

7.8 

CI2 

- 

4.0 

- 

N2 

6.7 

8.0 

- 

H2O 

9.3 

8.8 

- 

a.  Rate  constants  in  units  of  10-iO  cm3  s-1.  The  estimated  uncertainties  are 
±0.5x10-10  cm3  s-1. 


dependence  were  identified.  Spectra  used  to  characterize  energy  transfer  were  then 
recorded  at  pressure-delay  combinations  that  were  well  below  those  that  revealed  the  effects 
of  multiple  collisions. 

B.  Data  analysis 

Computer  simulations  of  the  pump-dump-probe  spectra  were  used  to  estimate  the 
state-to-state  energy  transfer  rate  constants.  The  computer  model  was  based  on  a  kinetic 
analysis  that  assumed  single  collision  conditions  and  that  collisions  with  a  single  species, 
M,  dominated.  Starting  with  all  of  the  molecules  in  the  prepared  level  vi  Ji,  the  equations 
that  describe  the  time  evolution  of  the  ground  state  population  distribution  are 

=ki[viJi],  and  =kf[viJi] 

where  [vi  JJ  and  [vf  Jf]  are  the  populations  in  the  initial  and  final  states,  respectively. 

Constants  ki  and  kf  are  defined  by  the  equations 

ki  =  I  IkM(viJi->vf  Jf)[M]  and  kf  =  kwCvi  Ji ->  vf  Jf)[M] 

Vf  Jf 

where  kM(vi  Ji  ->  vf  Jf)  is  a  state-to-state  transfer  rate  constant.  Integration  of  the  rate 
equations  yields  the  results 

[Vi  Ji]  =  [Vi  Jilo  e-kit  and  [Vf  Jf]  =  [vi  JJq  |  d  '  e’kit) 

Ro-vibrational  population  distributions,  calculated  from  the  above  equations  and 
trial  values  for  the  rate  constants,  were  used  to  generate  synthetic  spectra  that  could  be 
compared  with  the  experimental  results.  In  principle,  all  of  the  rate  constants  could  be 
varied  until  an  optimum  fit  was  achieved.  However,  the  quality  of  our  data  was  not 
sufficient  to  warrant  this  approach.  Consequently,  we  used  scaling  laws  to  parameterize 
the  rate  constants,  and  then  adjusted  the  parameter  values  to  obtain  agreement  with 
experiment.  For  analysis  of  the  pure  rotational  transfer  data,  a  simplified  version  of  the 
statistical  power  gap  law^^  (SPG)  was  used,  vis 

kM(Ji-^Jf)  =  CN^IAE/BrT^ 

(in  this  equation  the  vi=vf  labels  have  been  omitted).  C  and  y  are  fitted  parameters.  B  is 

the  rotational  constant,  AE  is  the  difference  in  energy  between  the  initial  and  final  states, 
and  is  a  statistical  factor  that  depends  on  the  extent  to  which  the  angular  momentum 

projection  quantum  number,  m,  is  modified  by  collisions.  On  the  basis  of  on  this 
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relationship,  graphical  simulations  were  used  to  obtain  best  estimates  for  the  parameters  C 

and  y.  The  C  coefficients  were  constrained  to  give  total  rotational  transfer  rate  constants 

(k^j  =  X  kM(Ji  ->  Jf))  that  were  compatible  with  the  decay  rate  constants  given  in  Table  3. 
Jf 

One  set  of  fittings  was  performed  with  the  assumption  that  m  is  randomized  by  collision.  A 
second  set  was  made  with  the  constraint  that  Am=0.  We  found  that  transfer  induced  by  the 
atomic  collision  partners  was  consistent  with  Am=0.  For  the  molecular  collision  partners, 
fits  obtained  with  the  assumption  of  m-randomizing  events  were  slightly  better  than  those 
obtained  with  Am=0.  Rotational  energy  transfer  parameters  are  collected  in  Table  4. 

Vibrational  energy  transfer  was  typically  an  order  of  magnitude  slower  than 
rotational  energy  transfer.  In  the  studies  conducted  to  date,  we  have  only  observed  Av=-1 
events.  Upward  transfer  or  multi-quantum  downward  transfer  was  not  detected.  Rotational 
contours  of  the  bands  originating  from  the  Av=-1  levels  reflected  a  tendency  to  conserve 
angular  momentum.  From  trial  simulations  it  was  found  that  rate  constants  defined  by  the 
expression 

knCvi  Ji— >Vi-l  Jf)  =  Cy^v-l  Na  lAErot/BI  ^ 

yielded  distributions  that  were  in  good  agreement  with  experiment.  Here,  AErot  is  the 
difference  in  rotational  energy  between  the  initial  and  final  states.  Note  that  the  equation  has 
a  singularity  at  AErot=0.  This  problem  was  avoided  by  setting  kwCvi  Ji->vi.i  Ji)  =  kM(vi  Ji 
^vi-i  Ji+2).  The  values  of  y  used  in  these  simulations  were  taken  from  Table  4,  leaving 
only  the  Cv-»v-l  coefficients  to  be  determined  from  the  fits.  Overall  rate  constants  for 
transfer  from  vi  to  vi-1  were  evaluated  by  summing  the  state-to-state  rate  constants. 
Vibrational  transfer  rate  constants  are  given  in  Table  5. 

C.  Modeling  ofl2^  deactivation 

After  the  slow  initiation  stage,  the  primary  channel  for  excitation  of  l2^  in  COIL 
systems  is  I2(X)  +  I(^Pi/2)  — >  l2^  +  I(^P3/2)-  We  have  shown  that  levels  in  the  range  of 
25<v"<46  are  populated  by  this  process.  The  energy  carried  by  I(^Pi/2)  or  02(^A) 
metastables  is  sufficient  to  dissociate  I2  from  levels  with  v">20.  However,  David  et  al.9 
have  suggested  that  only  levels  v"=30  and  higher  are  efficiently  dissociated,  due  to  Franck- 
Condon  restrictions.  In  modeling  the  deactivation  process,  we  have  treated  the  level  that 
defines  the  threshold  for  efficient  dissociation  (vth)  as  a  variable  parameter.  Hence,  the  rate 
at  which  l2^  is  deactivated  corresponds  to  the  rate  at  which  the  nascent  distribution  (c.f. 
Figure  4)  is  relaxed  to  levels  below  vth- 
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Rotational  Energy  Transfer  Rate  Constants  and  Statistical  Power  Gap  Law  Model  Parameters  for  I2(X). 


ON 


00 

cn 


m 


cn 

CN 


X) 

< 


ca 

O 


X 

< 

Z 


o 


t. 

t 

t 

CN 

CN 

CN 

CN 

’  i 

CN 

CN 

7t- 

m 

VO 

VO 

m 

VO 

VO 

d 

d 

d 

d 

d 

•rf 

Os 


Tl-  ^  O  ’-J  O') 

rt  in  so  vd  cn  NO 


t  t  t 

c3  cl  rl 


o  cn 

VO  VO 

d  d  d 


CN  t-; 

in  in  0C3  I  I  I 


I  - 

^  t:  ffi  < 

U  D- 


CS 


o 

^  ,  cs 

O  U  ffi 


IN  ^ 


o 

c« 

13 

> 

13 

c  • 

O  ’-H 

i  - 

£  S 


cn 

c 

o 

"c« 


O 

o 

00 

c 

‘E 

(U 

«3 

c 

o 

o 

I 

6 

XI 

'E 

73 

a> 

-2 

'o 

c 

c/s 

cd 


o 

.cd 


cd 

u 

CO 

cd 

CO 


o 

CO 

o 

•  w^ 

O 

1 

o 

C4-H 

Cm 

•— j 

o 

O 

Im 

o 

CO 

Uh 

^-H 

3 

Cm 

CO 

3 

O 

s 

3 

o 

W4 

p 

3 

a 

Im 

a> 

> 

O 

> 

4h 

•  w^ 

M 

00 

2a 

0) 

X 

g 

CO 

4*^ 

CO 

CO 

cn 

6 

3 

O 

1 

o 

X 

c.> 

•M 

o 

<4-4 

t/3 

g 

Cd 

o 

V/ 

U 

V 

cd 


Table  5 


Vibrational  Energy  Transfer  Rate  Constants  (Ad=-1)  for  l2(X). 


CoUision 

partner 

k(n"=23)a 

k(n"=38)a 

It 

1.2 

1.0 

Ar 

0.3 

0.4 

l2 

0.5 

0.6 

02 

” 

0.7 

Cl2 

— 

0.9 

H2O 

1.7 

2.6 

D2O 

— 

2.1 

a.  Rate  constants  in  units  of  lO’lO  cm^  s'l. 


The  relationship  between  the  measured  state-to-state  vibrational  energy  transfer  rate 
constants  and  the  deactivation  rate  constant  has  been  examined  by  computer  modeling.  The 
model  was  based  on  the  following  assumptions: 

(i)  Vibrational  energy  transfer  is  controlled  by  the  vf-vi=±l  selection  rule. 

(ii)  The  rate  constants  for  downward  transfer  ( )  are  dependent  on  the  vibrational 

quantum  number  in  a  simple  manner,  e.g.  linear  (Landau-Teller)  scaling;  =  vkj_^(,, 

or  power  law  scaling;  =  v^kj^o . 

(iii)  Upward  energy  transfer  rate  constants  can  be  calculated  from  detailed  balance; 

exp(-(E^  -  E^_,)  /  k  J). 

(iv)  The  rotational  distributions  are  thermalized  and  do  not  need  to  be  considered. 

The  initial  (time=0)  distribution  of  population  among  the  l2(X)  vibrational  levels 
was  defined  using  the  data  from  Figure  4.  The  number  densities  were  arbitrarily  chosen  to 
give  a  total  number  density  (summed  over  all  vibrational  levels)  of  3x10^0  cm'^  (the 
relaxation  kinetics  do  not  depend  on  the  I2  pressure  in  this  model). 

Coupled  rate  equations  that  describe  the  vibrational  energy  transfer  processes  are 
defined  by 

=  (k,_,^ai2,v  - 1]  +  k,,,_ JIj,v 4- 1] - -f-  k,^,_,)[l2, v])[M] 

These  equations  were  solved  by  numerical  integration.  At  equally  spaced  time  intervals 

59 

the  number  density  was  obtained  from  the  summation  [  ]=  E[I:,v). 

V=V|„ 

Figure  6  shows  a  typical  [  Ij]  vs  time  profile  for  relaxation  by  He.  This  calculation 
was  performed  with  linear  v-scaling  (n=l),  ko«_i=3.9xlO'^2  pnaS  g-l  (extrapolated  from  the 
SEP  data),  an  He  pressure  of  5  Torr,  and  T=298  K.  As  expected,  the  [l^]  time 
dependence  was  a  "sigmoid-like"  function;  [I^]  does  not  begin  to  decay  noticeably  until  a 
significant  fraction  of  the  population  has  relaxed  to  levels  near  v=20.  Characteristic 
lifetimes  for  these  curves  were  defined  by 

Mlctdt 
t  =  -2 - 

oe 

jUllc  dt 
0 

where  [I2]c=[l2]  ‘  fi2k=oo,  and  the  thermal  population  of  the  v>vth  levels  at 

equilibrium.  For  the  curve  shown  in  Figure  6,  this  expression  gives  T=  1.2  |is. 

To  begin  with,  the  model  was  used  to  examine  the  possibility  that  the  overall 
deactivation  process  could  be  represented  using  a  single  effective  rate  constant  (k^)  for 


10 


each  collision  partner.  Plots  of  1/x  vs.  collider-gas  pressure  were  found  to  be  linear, 
thereby  demonstrating  the  validity  of  the  concept  of  a  deactivation  rate  constant. 

Deactivation  rate  constants  were  calculated  from  the  measured  vibrational  energy 
transfer  rate  constants  for  He,  Ar,  CI2, 02,  and  H2O.  The  results  are  collected  in  Table  6. 
Deactivation  rate  constants  are  given  for  the  limiting  cases  of  vth=20  and  vth=30.  The 
power  scaling  (n)  used  for  each  calculation  was  chosen  in  accordance  with  the  available 
data.  For  He,  the  measurements  suggest  that  the  v->v-l  rate  constants  are  almost 
independent  of  v  over  the  range  20<v"<40.  However,  Table  6  lists  both  n=0  and  n=l 
results  for  He,  in  order  to  illustrate  the  effect  on  of  increasing  n  from  0  to  1.  These 
calculations  show  that,  provided  the  base  rate  constant  (k]_^o)  is  consistent  with  the 

measured  rate  constants,  changing  n  from  0  to  1  does  not  greatly  influence  the  deactivation 
rate  constants. 

D.  Significance  of  Deactivation  Rate  Constants  for  COIL  Models. 

Previous  models  of  COIL  systems  have  used  a  standard  set  of  estimated  rate 
constants  (AFWL  rate  package)  1^.  The  present  results  require  one  significant  change  in 
this  data  base.  The  AFWL  package  uses  a  rate  constant  for  deactivation  of  by  H2O  of 

3xl0'lb  cm3  s'L  Based  on  the  results  given  in  Table  6,  we  recommend  a  value  of 
kHp=1.5xlO'll  cm3  g-l 

The  rate  constant  for  deactivation  by  O2  was  estimated  at  5x10'^  1  cm3  5-1^  and  the 
present  results  are  consistent  with  this  value. 

As  He  is  used  as  the  carrier  gas  in  most  high-power  COIL  devices,  the  deactivation 
rate  constant  for  He  is  a  critical  parameter.  The  AFWL  rate  package  uses  kH^=4xl0'12  cm3 

s‘l,  while  the  present  results  indicate  a  slightly  higher  value,  within  the  range 
0.5<kH^<1.0xl0'^l  cm3  g-l.  Carrol^^  and  co-workers  have  suggested  that  kjj^  may  be 

larger  still.  Beyond  an  optimal  pressure,  COIL  output  power  decreases  with  increasing 
pressure.  Carrol^^  examined  the  possibility  that  deactivation  by  He  was  primarily 

responsible  for  this  behavior.  In  order  to  duplicate  the  observed  pressure  dependence, 
CarroH^  found  that  k^^  needed  to  be  increased  to  4.8xl0‘ll  cm3  when  this  was  the 

only  rate  constant  allowed  to  deviate  from  the  AFWL  values.  These  calculations  were 
performed  before  data  from  the  energy  transfer  measurements  were  available.  Since  the 
fitted  value  for  k^^  was  not  consistent  with  the  present  results,  a  model  with  two  variable 
rate  constants  is  being  investigated.  In  the  latest  calculations,  k^^  and  the  rate  constant 
(k25)  for  -I-  02(a^A)  ^  I  + 1  +  02(3E)  are  allowed  to  vary,  and  a  corrected  value  for 
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Table  6 


Effective  Deactivation  Rate  Constants  (298  K) 


Collision 

partner 

n 

ki->o^ 

He 

0 

12 

He 

1 

0.39 

Ar 

1 

0.12 

02 

0 

7.0 

Cl2 

0 

9.0 

H2O 

1 

0.71 

kt(vth=20)b  kt(vth=30)b 


5.4 

8.3 

5.1 

9.8 

1.6 

2.9 

3.3 

4.9 

4.1 

6.3 

9.3 
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a.  Units  of  10"11  cm3  s-1. 


b.  Units  of  10"12  cm3  g-l. 


is  used.  Respectable  fits  to  the  pressure  dependence  have  been  achieved^^  using 
k25=1.8xlO'10  cm3  g-l  and  kHg=1.6xlO"ll  cm3  g-l,  jhe  former  is  reduced  from  the 

AFWL  value  by  a  factor  of  0.6.  In  the  absence  of  experimental  data  for  k25,  this 
adjustment  seems  reasonable.  By  reducing  k25,  k^^  is  brought  much  closer  to  the  energy 

transfer  estimate.  It  is  probable  that  even  better  agreement  could  be  attained  by  treating 
more  of  the  uncertain  rate  constants  as  variables. 

Deactivation  of  I2  by  CI2  was  not  considered  in  the  AFWL  package  as  it  was 

thought  that  CI2  was  completely  consumed  in  the  02(a^A)  generator.  Recent  experiments 
suggest  that  Ci2  is  not  fully  utilized  when  02(alA)  generators  are  run  at  high  pressures^S. 
Furthermore,  there  was  speculation  that  rapid  deactivation  of  by  CI2  could  be  one  of  the 

factors  that  caused  the  loss  of  power  seen  at  high  pressures.  The  energy  transfer 
measurements  show  that  CI2  is  not  particularly  effective  in  deactivating  .  If  the  presence 

of  CI2  is  causing  a  loss  of  laser  power,  it  does  not  appear  to  operate  by  suppressing  the 
dissociation  of  I2. 

3.  Spectroscopy  of  Iodine  Monobromide. 

The  A  l(3n)  and  A'  2(3n)  states  of  IBr  were  observed  in  a  free-jet  expansion  of 

IBr  in  Ar,  following  photolysis  with  ArF  laser  irradiation. 

Laser  excitation  spectra  were  recorded  for  the  D'  2(3P2)  4-  A'  transition.  The  data 
included  high-resolution,  rotationaUy  analyzed  measurements  for  34  bands  of  and 

37  bands  of  plus  bandhead  measurements  for  >300  bands  recorded  in  low- 

resolution.  The  data  span  v"  =  0-5  and  v'  =  24-97.  Results  from  the  analysis  of  these 
data  were  merged  with  preliminary  results  from  a  study  of  the  D'  -»  A'  emission  spectrum 
to  yield  spectroscopic  constants  valid  for  v'  =  0-97,  v"  =  0-22.  The  primary 
spectroscopic  constants  (cm'l)  were:  ATg  =  28  018.98,  cOg'  =  123.088,  CDgXg'  =  0.2820, 
Bg’  =  0.02853,  cOg"  =  147.633,  cOgXg”  =  1.258,  B/'  =  0.04298.  This  work  is  described 
in  reference  29. 

Excitation  spectra  were  obtained  for  the  P  1(3112)  <—  A  transition  in  medium  (0.5 
cm‘1)  and  high  (0.08  cm'l)  resolution.  Sixteen  'u'-'o"  bands  recorded  in  high  resolution 
were  rotationaUy  analyzed,  providing  the  first  detailed  characterization  of  the  lowest  three  t) 
levels  in  the  A  state.  An  additional  50  bandheads  of  the  d"=0  and  1  progressions  were 
recorded  for  both  and  I^lBr.  The  new  data,  which  sample  "U'  levels  in  the  31-71 
range,  were  combined  with  existing  information  for  low  D  in  the  P  state  to  provide 
improved  constants  for  the  P  state.  These  results  were  merged  with  least-squares 
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parameters  from  a  reanalysis  of  existing  data  for  the  A  and  X  states  to  yield  spectroscopic 
parameters  valid  for  d=0-20  in  the  X  state  and  'i)=0-26  in  the  A  state.  Among  these  results 

are  the  first  precise  estimates  of  the  equilibrium  parameters  for  the  A  state: 
Te=12369.68(14)  cm'l,  C0e=l 34.22(1 3)  cm-1,  Be=0.04240(5)  cm'l,  Re=2.8583  A.  This 
work  is  described  in  reference  30. 
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Appendix 


RKR  Turning  Points  for  12(D) 
and 

Franck-Condon  Factors  for  the  D-X  Transition 
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Table  A1 


Rotational  Constants,  Vibrational  Energies,  and  RKR  Turning  Points  for  12(D) 


V 

B/cm-1 

G(v)/cm*l 

Rmax/^ 

0. 

0.0206930 

47.5037 

3.509 

3.658 

1. 

0.0206491 

142.3323 

3.458 

3.717 

2. 

0.0206050 

236.9427 

3.424 

3.759 

3. 

0.0205607 

331.3305 

3.397 

3.794 

4. 

0.0205163 

425.4916 

3.375 

3.825 

5. 

0.0204717 

519.4220 

3.355 

3.853 

6. 

0.0204270 

613.1179 

3.337 

3.880 

7. 

0.0203821 

706.5756 

3.321 

3.904 

8. 

0.0203370 

799.7916 

3.306 

3.928 

9. 

0.0202918 

892.7626 

3.292 

3.950 

10. 

0.0202465 

985.4853 

3.279 

3.971 

11. 

0.0202010 

1077.9566 

3.267 

3.992 

12. 

0.0201554 

1170.1737 

3.255 

4.012 

13. 

0.0201096 

1262.1336 

3.244 

4.032 

14. 

0.0200638 

1353.8337 

3.233 

4.051 

15. 

0.0200178 

1445.2714 

3.223 

4.070 

16. 

0.0199717 

1536.4443 

3.214 

4.088 

17. 

0.0199254 

1627.3500 

3.205 

4.106 

18. 

0.0198791 

1717.9863 

3.196 

4.124 

19. 

0.0198327 

1808.3511 

3.187 

4.141 

20. 

0.0197861 

1898.4423 

3.179 

4.158 

21. 

0.0197395 

1988.2580 

3.171 

4.175 

22. 

0.0196927 

2077.7963 

3.163 

4.192 

23. 

0.0196459 

2167.0556 

3.156 

4.208 

24. 

0.0195990 

2256.0341 

3.149 

4.225 

25. 

0.0195520 

2344.7303 

3.142 

4.241 

26. 

0.0195049 

2433.1427 

3.135 

4.257 

27. 

0.0194577 

2521.2700 

3.128 

4.272 

28. 

0.0194105 

2609.1107 

3.121 

4.288 

29. 

0.0193632 

2696.6636 

3.115 

4.304 

30. 

0.0193159 

2783.9275 

3.109 

4.319 

31. 

0.0192685 

2870.9014 

3.103 

4.334 

32. 

0.0192210 

2957.5841 

3.097 

4.350 

33. 

0.0191735 

3043.9747 

3.091 

4.365 

34. 

0.0191259 

3130.0722 

3.086 

4.380 

35. 

0.0190783 

3215.8759 

3.080 

4.395 

36. 

0.0190306 

3301.3848 

3.075 

4.409 

37. 

0.0189829 

3386.5983 

3.069 

4.424 

38. 

0.0189352 

3471.5156 

3.064 

4.439 

39. 

0.0188875 

3556.1361 

3.059 

4.454 

40. 

0.0188397 

3640.4593 

3.054 

4.468 
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Table  A1  -  continued 


V 

B/cm‘* 

G(v)/cm'^ 

Rmin/^ 

41. 

0.0187919 

3724.4845 

3.049 

42. 

0.0187441 

3808.2114 

3.044 

43. 

0.0186963 

3891.6394 

3.040 

44. 

0.0186484 

3974.7682 

3.035 

45. 

0.0186006 

4057.5973 

3.031 

46. 

0.0185528 

4140.1265 

3.026 

47. 

0.0185049 

4222.3556 

3.022 

48. 

0.0184571 

4304.2842 

3.017 

49. 

0.0184093 

4385.9121 

3.013 

50. 

0.0183615 

4467.2393 

3.009 

51. 

0.0183137 

4548.2656 

3.005 

52. 

0.0182660 

4628.9909 

3.001 

53. 

0.0182183 

4709.4152 

2.997 

54. 

0.0181706 

4789.5383 

2.993 

55. 

0.0181229 

4869.3604 

2.989 

56. 

0.0180753 

4948.8815 

2.985 

57. 

0.0180277 

5028.1017 

2.981 

58. 

0.0179802 

5107.0209 

2.978 

59. 

0.0179327 

5185.6395 

2.974 

60. 

0.0178853 

5263.9575 

2.970 

61. 

0.0178379 

5341.9751 

2.967 

62. 

0.0177906 

5419.6925 

2.963 

63. 

0.0177434 

5497.1100 

2.960 

64. 

0.0176962 

5574.2279 

2.956 

65. 

0.0176491 

5651.0463 

2.953 

66. 

0.0176021 

5727.5657 

2.949 

67. 

0.0175552 

5803.7863 

2.946 

68. 

0.0175084 

5879.7085 

2.943 

69. 

0.0174616 

5955.3327 

2.940 

70. 

0.0174150 

6030.6592 

2.936 

71. 

0.0173684 

6105.6886 

2.933 

72. 

0.0173220 

6180.4212 

2.930 

73. 

0.0172756 

6254.8575 

2.927 

74. 

0.0172294 

6328.9980 

2.924 

75. 

0.0171833 

6402.8431 

2.921 

76. 

0.0171373 

6476.3934 

2.918 

77. 

0.0170914 

6549.6495 

2.915 

78. 

0.0170456 

6622.6118 

2.912 

79. 

0.0170000 

6695.2810 

2.909 

80. 

0.0169545 

6767.6577 

2.906 

81. 

0.0169092 

6839.7424 

2.903 

82. 

0.0168639 

6911.5357 

2.900 

83. 

0.0168189 

6983.0384 

2.897 

Rmax/^ 

4.483 

4.497 

4.512 

4.526 

4.540 

4.554 

4.569 

4.583 

4.597 

4.611 

4.625 

4.639 

4.653 

4.667 

4.681 

4.695 

4.709 

4.723 

4.737 

4.751 

4.765 

4.779 

4.793 

4.806 

4.820 

4.834 

4.848 

4.862 

4.875 

4.889 

4.903 

4.917 

4.930 

4.944 

4.958 

4.972 

4.986 

4.999 

5.013 

5.027 

5.041 

5.054 

5.068 
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Table  A1  -  continued 


V 

B/cm-1 

G(v)/cm'^ 

R-min/^ 

Rmax/^ 

84. 

0.0167740 

7054.2511 

2.895 

5.082 

85. 

0.0167292 

7125.1743 

2.892 

5.096 

86. 

0.0166846 

7195.8089 

2.889 

5.109 

87. 

0.0166402 

7266.1556 

2.886 

5.123 

88. 

0.0165959 

7336.2149 

2.884 

5.137 

89. 

0.0165518 

7405.9878 

2.881 

5.151 

90. 

0.0165078 

7475.4748 

2.878 

5.165 

91. 

0.0164641 

7544.6769 

2.876 

5.178 

92. 

0.0164205 

7613.5947 

2.873 

5.192 

93. 

0.0163771 

7682.2290 

2.871 

5.206 

94. 

0.0163339 

7750.5806 

2.868 

5.220 

95. 

0.0162909 

7818.6504 

2.865 

5.234 

96. 

0.0162481 

7886.4392 

2.863 

5.247 

97. 

0.0162055 

7953.9478 

2.860 

5.261 

98. 

0.0161631 

8021.1770 

2.858 

5.275 

99. 

0.0161209 

8088.1278 

2.855 

5.289 

100. 

0.0160789 

8154.8009 

2.853 

5.303 

101. 

0.0160372 

8221.1973 

2.850 

5.317 

102. 

0.0159956 

8287.3179 

2.848 

5.330 

103. 

0.0159543 

8353.1635 

2.845 

5.344 

104. 

0.0159132 

8418.7352 

2.843 

5.358 

105. 

0.0158724 

8484.0337 

2.840 

5.372 

106. 

0.0158318 

8549.0601 

2.838 

5.386 

107. 

0.0157914 

8613.8153 

2.836 

5.400 

108. 

0.0157512 

8678.3003 

2.833 

5.414 

109. 

0.0157114 

8742.5159 

2.831 

5.428 

110. 

0.0156717 

8806.4633 

2.828 

5.442 

111. 

0.0156323 

8870.1433 

2.826 

5.455 

112. 

0.0155932 

8933.5570 

2.824 

5.469 

113. 

0.0155544 

8996.7054 

2.821 

5.483 

114. 

0.0155158 

9059.5894 

2.819 

5.497 

115. 

0.0154775 

9122.2101 

2.817 

5.511 

116. 

0.0154395 

9184.5685 

2.814 

5.525 

117. 

0.0154017 

9246.6657 

2.812 

5.539 

118. 

0.0153642 

9308.5026 

2.810 

5.553 

119. 

0.0153270 

9370.0803 

2.807 

5.567 

120. 

0.0152902 

9431.3999 

2.805 

5.581 

121. 

0.0152536 

9492.4625 

2.803 

5.595 

122. 

0.0152173 

9553.2690 

2.800 

5.609 

123. 

0.0151813 

9613.8205 

2.798 

5.623 

124. 

0.0151456 

9674.1182 

2.796 

5.637 

125. 

0.0151102 

9734.1631 

2.794 

5.651 

20 


Table  A2 

Franck-Condon  Factors  for  the  D-X  System  of  I2 


v" 

v'=  4 

II 

v’=  6 

II 

21 

8.0573E-14 

6.2078E-13 

4.0116E-12 

2.2360E-11 

22 

5.5414E-13 

4.1028E-12 

2.5471E-11 

1.3634E-10 

23 

3.5730E-12 

2.5380E-11 

1.5 11  IE- 10 

7.7540E-10 

24 

2.1620E-11 

1.4707E-10 

8.3819E-10 

4.1152E-09 

25 

1.2285E-10 

7.9879E-10 

4.3490E-09 

2.0385E-08 

26 

6.5592E-10 

4.0677E-09 

2.1109E-08 

9.4246E-08 

27 

3.2912E-09 

1.9422E-08 

9.5835E-08 

4.0649E-07 

28 

1.5521E-08 

8.6932E-08 

4.0676E-07 

1.6343E-06 

29 

6.8774E-08 

3.6459E-07 

1.6128E-06 

6.1185E-06 

30 

2.8624E-07 

1.4317E-06 

5.9674E-06 

2.1297E-05 

31 

1.1182E-06 

5.2586E-06 

2.0572E-05 

6.878  lE-05 

32 

4.0958E-06 

1.8038E-05 

6.5950E-05 

2.0557E-04 

33 

1.4048E-05 

5.768  lE-05 

1.9610E-04 

5.6667E-04 

34 

4.5045E-05 

1.7154E-04 

5.3915E-04 

1.4347E-03 

35 

1.3473E-04 

4.7302E-04 

1.3650E-03 

3.3179E-03 

36 

3.7486E-04 

1.2049E-03 

3.1659E-03 

6.9592E-03 

37 

9.6703E-04 

2.8217E-03 

6.6822E-03 

1.3115E-02 

38 

2.303  lE-03 

6.0383E-03 

1.2723E-02 

2.1923E-02 

39 

5.0369E-03 

1.1714E-02 

2.1590E-02 

3.1903E-02 

40 

1.0045E-02 

2.0379E-02 

3.2106E-02 

3.9282E-02 

41 

1.8100E-02 

3.1326E-02 

4.0777E-02 

3.8988E-02 

42 

2.9103E-02 

4.1618E-02 

4.2394E-02 

2.8314E-02 

43 

4.1009E-02 

4.6133E-02 

3.3264E-02 

1.1591E-02 

44 

4.9270E-02 

4.0035E-02 

1.6119E-02 

4.7852E-04 

45 

4.8172E-02 

2.3607E-02 

1.8991E-03 

5.0983E-03 

46 

3.4945E-02 

5.8136E-03 

2.7669E-03 

2.1071E-02 

47 

1.4715E-02 

4.3360E-04 

1.8517E-02 

2.927  lE-02 

v" 

v’=  8 

v'=  9 

v'=  10 

v’=  11 

21 

1.0972E-10 

4.8145E-10 

1.9122E-09 

6.9424E-09 

22 

6.4228E-10 

2.7043E-09 

1.0302E-08 

3.5859E-08 

23 

3.500  lE-09 

1.4114E-08 

5.1468E-08 

1.7137E-07 

24 

1.7763E-08 

6.8455E-08 

2.3840E-07 

7.5760E-07 

25 

8.3952E-08 

3.0847E-07 

1.0234E-06 

3.0955E-06 

26 

3.6940E-07 

1.2907E-06 

4.0679E-06 

1.1676E-05 

27 

1.5122E-06 

5.009  lE-06 

1.4950E-05 

4.0580E-05 

28 

5.7524E-06 

1.8006E-05 

5.0704E-05 

1.2963E-04 

29 

2.0305E-05 

5.9826E-05 

1.5829E-04 

3.7942E-04 

30 

6.6370E-05 

1.8326E-04 

4.5332E-04 

1.0132E-03 

31 

2.0036E-04 

5.1578E-04 

1.1859E-03 

2.4547E-03 

32 

5.5670E-04 

1.3279E-03 

2.8177E-03 

5.3568E-03 

33 

1.4175E-03 

3.1095E-03 

6.0357E-03 

1.0427E-02 

21 


Table  A2  -  Continued 


v" 

v’=  8 

v'=  9 

v'=  10 

v'=  11 

34 

3.289  lE-03 

6.5740E-03 

1.1541E-02 

1.7860E-02 

35 

6.9040E-03 

1.2423E-02 

1.9427E-02 

2.6396E-02 

36 

1.2981E-02 

2.0697E-02 

2.8217E-02 

3.2647E-02 

37 

2.1570E-02 

2.9790E-02 

3.4267E-02 

3.2053E-02 

38 

3.1059E-02 

3.5901E-02 

3.2945E-02 

2.2419E-02 

39 

3.7589E-02 

3.4295E-02 

2.2381E-02 

8.2013E-03 

40 

3.6273E-02 

2.3171E-02 

7.6913E-03 

8.0793E-05 

41 

2.5049E-02 

7.9139E-03 

1.5700E-05 

5.7884E-03 

42 

9.0798E-03 

1.3433E-05 

6.7283E-03 

1.8901E-02 

43 

8.4013E-05 

7.0523E-03 

2.0366E-02 

2.2252E-02 

44 

6.5588E-03 

2.1386E-02 

2.3156E-02 

9.9762E-03 

45 

2.1771E-02 

2.4604E-02 

9.961  lE-03 

2.1600E-05 

46 

2.6701E-02 

1.0946E-02 

2.8030E-06 

8.9155E-03 

47 

1.3301E-02 

2.4549E-05 

9.8185E-03 

2.0836E-02 

v" 

v’=  12 

v’=  13 

v'=  14 

v'=  15 

21 

2.3228E-08 

7.2103E-08 

2.0884E-07 

5.6718E-07 

22 

1.1496E-07 

3.4173E-07 

9.4728E-07 

2.4604E-06 

23 

5.2525E-07 

1.4917E-06 

3.9472E-06 

9.7784E-06 

24 

2.2142E-06 

5.9913E-06 

1.5091E-05 

3.5547E-05 

25 

8.6030E-06 

2.2111E-05 

5.2837E-05 

1.1792E-04 

26 

3.0755E-05 

7.4819E-05 

1.6897E-04 

3.5579E-04 

27 

1.0093E-04 

2.3146E-04 

4.9179E-04 

9.7208E-04 

28 

3.0313E-04 

6.5208E-04 

1.2964E-03 

2.3908E-03 

29 

8.2978E-04 

1.6646E-03 

3.0761E-03 

5.2517E-03 

30 

2.0593E-03 

3.8249E-03 

6.5142E-03 

1.0193E-02 

31 

4.6009E-03 

7.8395E-03 

1.2169E-02 

1.7215E-02 

32 

9.1672E-03 

1.4156E-02 

1.9725E-02 

2.473  lE-02 

33 

1.6078E-02 

2.2122E-02 

2.7056E-02 

2.9140E-02 

34 

2.4357E-02 

2.9120E-02 

3.0145E-02 

2.6360E-02 

35 

3.0961E-02 

3.085  lE-02 

2.5259E-02 

1.5791E-02 

36 

3.1423E-02 

2.4079E-02 

1.3238E-02 

3.7357E-03 

37 

2.3058E-02 

1.1046E-02 

1.9954E-03 

4.0020E-04 

38 

9.3408E-03 

9.0181E-04 

1.4736E-03 

9.1106E-03 

39 

3.2084E-04 

2.9184E-03 

1.2086E-02 

1.8935E-02 

40 

4.4429E-03 

1.4759E-02 

2.0296E-02 

1.5687E-02 

41 

1.7035E-02 

2.1113E-02 

1.3855E-02 

3.2898E-03 

42 

2.1665E-02 

1.2117E-02 

1.5988E-03 

1.5255E-03 

43 

1.0756E-02 

6.1266E-04 

3.3890E-03 

1.3383E-02 

44 

1.6042E-04 

5.4562E-03 

1.5964E-02 

1.6026E-02 

45 

7.3862E-03 

1.7952E-02 

1.4900E-02 

3.5048E-03 

46 

1.9500E-02 

1.3661E-02 

1.7259E-03 

2.3336E-03 

47 

1.2706E-02 

6.8912E-04 

4.6028E-03 

1.5090E-02 

22 


Table  A2  -  Continued 


v" 

v'=  16 

v'  =  17 

v’=  18 

v'=  19 

21 

1.4503E-06 

3.5047E-06 

8.0287E-06 

1.7485E-05 

22 

6.0125E-06 

1.3873E-05 

3.0320E-05 

6.2930E-05 

23 

2.2770E-05 

5.0013E-05 

1.0393E-04 

2.0484E-04 

24 

7.8614E-05 

1.6378E-04 

3.2233E-04 

6.0072E-04 

25 

2.467  lE-04 

4.8542E-04 

9.0057E-04 

1.5787E-03 

26 

7.0095E-04 

1.2958E-03 

2.253  lE-03 

3.6907E-03 

27 

1.7932E-03 

3.0949E-03 

5.0063E-03 

7.5986E-03 

28 

4.1006E-03 

6.5534E-03 

9.7686E-03 

1.3581E-02 

29 

8.2993E-03 

1.2150E-02 

1.6470E-02 

2.0632E-02 

30 

1.4661E-02 

1.9365E-02 

2.3415E-02 

2.5762E-02 

31 

2.2150E-02 

2.5802E-02 

2.6972E-02 

2.4907E-02 

32 

2.7713E-02 

2.7405E-02 

2.3364E-02 

1.641  lE-02 

33 

2.7147E-02 

2.1133E-02 

1.2775E-02 

4.9293E-03 

34 

1.8509E-02 

9.2554E-03 

2.1827E-03 

4.2879E-05 

35 

6.1765E-03 

5.3990E-04 

1.0802E-03 

6.6294E-03 

36 

5.3369E-07 

3.2261E-03 

1.0402E-02 

1.6170E-02 

37 

6.0449E-03 

1.3933E-02 

1.7806E-02 

1.4742E-02 

38 

1.683  lE-02 

1.8080E-02 

1.1956E-02 

3.6992E-03 

39 

1.7239E-02 

8.7488E-03 

1.1949E-03 

9.1570E-04 

40 

5.7268E-03 

6.7484E-05 

3.3838E-03 

1.1113E-02 

41 

2.7868E-04 

6.7148E-03 

1.4254E-02 

1.4201E-02 

42 

1.0217E-02 

1.6083E-02 

1.2053E-02 

3.3334E-03 

43 

1.6573E-02 

9.0729E-03 

8.8020E-04 

1.8643E-03 

44 

6.0382E-03 

2.0878E-06 

5.1061E-03 

1.2910E-02 

45 

6.2412E-04 

8.8789E-03 

1.4663E-02 

9.4754E-03 

46 

1.2362E-02 

1.463  lE-02 

5.9151E-03 

2.8129E-06 

47 

1.3234E-02 

2.7761E-03 

9.6820E-04 

9.2594E-03 

23 


